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Abstract 

The hypothesis that phase transitions originate from some topological change of the critical level 
hypersurface of the potential energy receives direct evident support by our study of the Bishop- 
Peyrard model of DNA thermal denaturation. 
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In spite of the huge body of hterature on the subject of thermodynamic phase transitions, 
one has at least two reasons to not be completely satisfied with our present understanding 
of such prominent phenomenon: From a theoretical viewpoint, the loss of analyticity in 
the derivatives of the free energy P| and the consequent appearance of singularities in the 
thermodynamic observables are only macroscopic signals of the phase transition, while the 
physical origin should be looked for at the microscopical dynamics level Q]. On the other 
hand, transitional phenomena experimented in finite atomic clusters, spin glasses, and bi- 
ological systems still challenge physicists to find out new tools and methods to cope with 
their difficulties. 

Both aspects of the problem are addressed by a recently proposed geometric approach to 
phase transitions 0,0,0]. Given any Hamiltonian system, the configuration space can be 
endowed with a metric, in order to obtain a Riemannian geometrization of the dynamics 0|. 
At the beginning, several numerical and analytical studies of a variety of models showed that 
the fiuctuation of the curvature becomes singular at the transition point. Then in Ref. 
the following conjecture was proposed: The phase transition is determined by a change in 
the topology of the configuration space, and the loss of analyticity in the thermodynamic 
observables is nothing but a consequence of such topological change. The latter conjecture 
is also known as the topological hypothesis. 

The aim of this Letter is to shed some other light upon the actual mechanism underlying 
the topological change. In order to do this, we choose to study a system which, though 
one dimensional (ID), has the peculiarity to exhibit a second order phase transition, the 
Peyrard-Bishop model ^l, developed as a simple model for describing the DNA thermally 
induced denaturation. Its simplicity already permitted an analytical computation of the 
largest Lyapunov exponent ^ by exploiting the geometrical method proposed in Ref. ^. 
On the other hand, we pursued a study of the topological characterization of the phase 
transition in a full analytic way. The Letter is organized as follows: First, we give a brief 
description of the main features of the Peyrard-Bishop model, the main body is devoted to 
the study of the topology of the configuration space of the model, and, last, we discuss our 
results and state our conclusions. 
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THE MODEL 

The system we studied was introduced in 1989 by Peyrard and Bishop '6] as a simple 
model for DNA thermal denaturation. It is defined by the following Hamiltonian 



N 



^2 



j=i 

which represents the energy of a string of base pairs of reduced mass m. Each hydrogen 
bond is characterized by the stretching i/i and its conjugate momentum pi = m{dyi/ dt). 
The elastic transverse force between neighboring pairs is tuned by the constant while the 
energy D and the inverse length a determine, respectively, the plateau and the narrowness 
of the on-site potential well that mimics the interaction between bases in each pair. It is 
understood that K, D, and a are all positive parameters. The transverse, external stress 
h > is a computational tool useful in the evaluation of the susceptibility. Our interest in 
it lies in the fact that a phase transition can occur only when h = 0. We assume periodic 
boundary conditions. 

The transfer operator technique 0, 0] maps the problem of computing the classical parti- 
tion function into the easier task of evaluating the lowest energy eigenvalues of a "quantum" 
mechanical Morse oscillator (no real quantum mechanics is involved, since the temperature 
plays the role of h). One can then observe that, as the temperature increases, the number 
of levels belonging to the discrete spectrum decreases, until for some critical temperature 



Tc = 2y2KD / {akB) only the continuous spectrum survives. This passage from a localized 
ground state to an unnormalizable one corresponds to the second order phase transition 
of the statistical model. Various critical exponents can be analytically computed and all 
applicable scaling laws can be checked. Detailed calculations as well as references to original 
works can be found in Ref. 



TOPOLOGY OF CONFIGURATION SPACE 

According to the Morse theory it is possible to understand the topology of a given 
manifold by studying the regular critical points of some smooth function (called Morse 
function) defined on it. In our case, the manifold M. is the configuration space and the 
natural choice for the Morse function is the potential V{y). Hence, one is lead to define the 
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following family of submanifolds of A4: 

M, = {yeR'': V{y) < v} . (2) 

A full characterization of the topological properties of Ai^ generally requires one to find the 
critical points of V{y), which means solving the equations 

dV 

Moreover, one has to compute the indexes of all the critical points, that is the number of 
negative eigenvalues of the Hessian d'^V/ (dyidyj). Then the Euler characteristic xi-^v) can 
be computed by means of the formula 

N 

xiM,) = ^(-l)Vfc(-M.) , (4) 

fc=0 

where ^k{M.v) is the total number of critical points of V{y) on M.^ which have index 
k. If applied to any generic model, such a task turns out to be quite formidable, but 
the exceptional simplicity of the Peyrard-Bishop model Q makes it possible to carry on 
completely the topological analysis without invoking Eq. (jlj). 

For the potential in exam, Eq. (jH)) results in the nonlinear system 

- 2y, + = h- 2{e-^^y^ - e''^^^ , (5) 

where R = Da^/K is a dimensionless ratio. It is easy to verify that a particular solution is 
given by 



1 i + VT+2h . 

yi = — In ^ = l,...,iV. (6) 



The corresponding minimum of potential energy is 



2 2 

Actually, there is only one other solution. Here is the proof: Equation can also be 
interpreted as a finite difference second order equation for the evolution of the variable yi. 
Let us introduce a set of new variables xi = {yi — yi-i). We can fix an initial condition 
Vn = Vo = y so that this definition can be reversed into yi = y + J2]=i ^j- terms of these 
new variables, Eq. (0) is written 



R 

a 



h — 2 exp ^ — 2ay — 2a xj^ + 2 exp ^ — ay — a xj^ 



(7) 



while the periodicity condition reads simply ^^^iXj = 0. Now we prove that no solution 
is compatible with this requirement besides the trivial solution Xj = Vz = 1, . . . , A^. From 
the recursive relation Eq. ((7j), we deduce that Xj+i > if and only if 



exp y-ay " « j ^ ^ • (8) 

Let us suppose that for some k it results Xk > 0. If Xk > Xk-i, then a fortiori Xk+i > Xk 
because e~"^'= < 1. Now Xk+i > and Xk+i > Xk imply 0:^4.2 > Xk+i and so on: The 
series {xj} increases for j > k. On the other hand, if Xk < Xk^i, then e~°'^'^-^ < 1; hence, 
Xk-i < Xk-2 and so on. Either situation is incompatible with a periodic boundary condition. 
A completely analogous argument rules out the possibility that it could be < for some 
value of k. 

We have proven that no solution is possible for Eq. ^ unless it results Ui = Y \Ii. Now 
it is evident that, when /i = 0, the limiting point Y — > +00 is a critical point. Here it 
is useful to introduce new variables Wi = e~"^% so that the condition of vanishing gradient 
reads 

1 

^In = h - 2w^iwi - 1) , (9) 

from which we immediately recognize, beside the minimum Eq. a second solution Wi = 
Vi = 1, . . . , A^, present only when h = 0. This is an irregular critical point because its 
Hessian matrix is degenerate; hence, the potential can no longer be considered a good Morse 
function. Two remarks are in order: First, as already stated, we do not need Eq. Q because 
My is homotopous to a full disk; second, since the Euler characteristic is a homotopical 
invariant, it does not change at Vc for the submanifolds Aiy 

Notwithstanding this, we still have a way to complete our analysis. Let us consider the 
boundaries of M^, the hypersurfaces 

S„ = {yeR'' : V{y) = v} . (10) 

If we choose one direction y = {iji, . . . , y^) and study the asymptotic behavior of the function 
V{y) in that direction by taking the limit A — > 00 in V{Xy), we find that if there is at least 
one negative component yj < 0, then the limit is positive infinity because of the exponential 
term. If all of the components are positive, then the potential goes as 

Y,im+i - m) + XDha Y,m + ND, (11) 
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FIG. 1: Level hyper surf aces, N = 2. The black line is the critical hypersurface. 



so that the limit is always positive infinity unless all yj are equal and /i = 0. This implies 
that for V raising from V^in to Vc = ND it happens that S^, is a closed (N-l) hypersurface 
[x(S^) = or 2 depending on N being odd or even, respectively]. On the contrary, if v = Vc 
and /i = 0, in the direction singled out by the condition yi = y2 = ■ ■ ■ = yN, the level 
hypersurface S^,^ closes at positive infinity (see Fig. „ for a = 2 example), while for each 
V > Vc the level hypersurfaces fail to close [x(S^) = 1]. This topological change is evidently 
related to the phase transition, since its presence is due to the very same feature of the 
potential (its plateau) giving rise to the phase transition. 



CONCLUSIONS 



Since its original formulation j^, the topological hypothesis has received several indirect 
and direct confirmations. If we limit ourselves to the latter, we may cite the numerical 
study of the 0^ theory Ql and the analytical computations of the Euler characteristic of the 

13] and of the mean- field fc-trigonometric model jl4|. 



mean-field XY model 



The main difference between our work and previous results lies in the low dimensionality 
of the Peyrard-Bishop model. By choosing a ID system, we get an enormous simplification 
of the landscape of critical points of the potential V{y), resulting in a completely trans- 
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parent description of the relation between topology and phase transition. Indeed, the same 
mechanism responsible for the appearance of the phase transition is clearly seen to produce 
the topology change in configuration space. 

Moreover, the present work contributes to enlarge the record of cases (quite a useful job 
in a research field so recent and full of open questions such as this) in at least two respects: 
It is the first analytical study of the topology of a noncompact configuration manifold and, 
while previous results 3] identify the phase transition point in a discontinuity of the first 
derivative of xi-M.^), the present work is peculiar in the sense that the topological change 
is signaled by the Euler characteristic of the boundary submanifolds xi^v)- 

In conclusion, it should be remarked that, though the Peyrard-Bishop model is very 
simple, it represents a first step in the exploration of the complex and interesting world 
of biological transition phenomena. Topological techniques will hopefully provide valuable 
insight to such a challenging research topic. 
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